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M ethod and system for processing a sound field representation 

The present invention relates to a method and a system for processing 
a sound field representation in order to provide a representation of that sound 
field modified by means of at least spatial processing. 

Methods and systems currently exist which, starting from a sound field 
5 representation in the form of coefficients as a function of time and space, enable 
a modified representation of those coefficients to be provided by linear 
combinations corresponding to at least spatial processing of the sound field. 

However, the possible processing operations are limited and are 
predetermined in all cases. 
10 For example, the systems grouped under the term « ambisonic » use a 

representation of a sound environment in the form of an angular distribution of 
virtual sources around a point assimilated to the listening position. 

Those systems permit the performance of only specific predetermined 
modifications, such as an angular distortion on the azimuth, focusing in a 
15 privileged direction, or rotation relative to a given axis. 

In addition, those systems are limited as regards their precision. Thus, 
in existing systems, processing operations, such as focusing in any direction, can 
be applied only to representations limited to a level of precision corresponding to 
first-order mathematics, that is to say, of poor quality. 
20 It therefore appears that none of the systems of the prior art permits 

the performance of any modification of a sound field representation in order to 
obtain a modified representation, especially in the case of an initial representation 
of excellent quality, for example of high order. 

The object of the invention is to overcome that problem by providing a 
25 method and a system for processing a sound field representation, permitting the 
performance of a very wide variety of processing operations. 

The invention relates to a method for processing a sound field 
representation in order to provide a representation modified by the application of 
this processing corresponding to at least spatial processing of the sound field, 
30 comprising: 

- the establishment of an initial representation of the sound field in the 
form of a set of coefficients representative of that field in time and in the three 
spatial dimensions; and 



- the establishment of a representation modified by filtering 
combinations that correspond to the processing operation and that are applied to 
the coefficients of the initial representation, 

characterized in that it also comprises a step of defining the processing 
operation, which step comprises: 

- a step of establishing a set of directivity functions which is 
representative of the processing operation, in the form of a set of coefficients 
corresponding to the decomposition of those functions on a basis of spherical 
harmonics; and 

- a step of determining the filtering combinations on the basis of that 
decomposition for the implementation of a step of applying the filtering 
combinations so defined to the coefficients of the initial representation, in order to 
obtain the modified representation. 

According to other features of the invention: 

- the step of determining a set of directivity functions comprises a sub- 
step for the direct determination of all or part of the set of directivity functions in 
one or other of a plurality of predetermined formats; 

- the step of determining a set of directivity functions comprises a sub- 
step for the automatic interpretation of a processing instruction permitting the 
automatic provision of all or part of the set of directivity functions in one or other 
of a plurality of predetermined formats, on the basis of the processing instruction; 

- the set of directivity functions which is representative of the 
processing operation is expressed in accordance with one or other of a plurality 
of predetermined formats and the step of establishing a set of directivity functions 
comprises a sub-step of conversion, from the predetermined format, into a set of 
coefficients corresponding to the decomposition of those directivity functions on a 
basis of spherical harmonics; 

- the step of defining the processing operation also comprises a step of 
establishing a set of parameters for that processing operation; 

- the set of parameters comprises all or some of the following 
parameters: 

- a set of constants forming weighting coefficients and defined as a 
function of the nature of the processing operation; and 
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- a parameter specifying the limit order of operation of the method 
and corresponding to the desired mathematical precision limit; 

- the step of establishing a set of parameters comprises a sub-step for 
the direct determination of all or part of the set of parameters; 

5 - the step of establishing a set of parameters comprises a sub-step for 

the automatic interpretation of a processing instruction permitting the automatic 
provision of all or part of the set of parameters on the basis of that processing 
instruction; 

- the step of establishing a set of parameters comprises a sub-step of 
10 calculating all or part of the set of constants as a function of the nature of the 

processing operation; 

- the step of determining the filtering combinations comprises at least 
one sub-step of calculating filtering combinations by means of linear 
combinations of the coefficients corresponding to the spherical harmonic-based 

15 decomposition of the set of directivity functions which is representative of the 
processing operation; 

- the sub-step of calculation by means of linear combinations uses the 
set of constants as weighting coefficients; 

- the procedure is repeated a plurality of times and the filtering 
20 combinations resulting from the various calculation sub-steps are combined with 

one another to form the filtering combinations that are to be applied to the 
coefficients of the initial representation; 

- the result of the calculation sub-step forms directly the filtering 
combinations that are to be applied to the coefficients of the initial representation; 

25 - the processing is convolution processing; 

- the processing is windowing processing; 

- the processing comprises at least one distortion to be applied to the 
initial sound field representation, and all or some of the directivity functions 
determined during the step of establishing a set of directivity functions form a set 

30 of distortion pairs which is representative of the distortion, in the form of a set of 
coefficients corresponding to the decomposition of the distortion pairs on a basis 
of spherical harmonics; 

- the processing comprises a first processing operation, which is 
referred to as nominal, and the distortion(s), the calculation sub-step providing 
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filtering combinations corresponding to the nominal processing and the step of 
determining filtering combinations also comprising at least one sub-step of 
distortion by modifying the filtering combinations provided by the calculation sub- 
step, on the basis of the distortion pairs representative of the distortion(s), in 
5 order to provide the filtering combinations that are to be applied to the coefficients 
of the initial representation; 

- the set of parameters comprises a parameter representative of the distortion 
effort, that is to say, specifying the resemblance between the processing to be 
applied and the nominal processing; 

10 - the processing comprises solely distortion(s), the step of 

determining the filtering combinations comprising at least one distortion sub-step 
providing filtering combinations solely on the basis of the distortion pairs 
representative of the distortion(s), in order to provide the filtering combinations 
that are to be applied to the coefficients of the initial representation; 

15 - it comprises a plurality of distortion sub-steps which are recursive 

with respect to one another, so that each distortion sub-step receives as an input 
the filtering combinations provided by the preceding sub-step; 

- the at least one distortion sub-step comprises an optimization sub- 
step permitting the determination of the filtering combinations corresponding to 

20 modifications to be carried out on portions of the initial representation that are not 
directly targeted by the distortion pairs; 

- the processing is split into a plurality of elemental processing 
operations and at least the step of determining filtering combinations is repeated 
a plurality of times, the filtering combinations resulting from the various 

25 calculation and/or distortion sub-steps and corresponding to the elemental 
processing operations being arranged amongst one another to form the filtering 
combinations that are to be applied to the coefficients of the initial representation; 

- the step of determining the filtering combinations comprises a sub- 
step of converting the filtering combinations represented in a predetermined 

30 format to provide filter parameters suitable for the parameterization of given filters 
with a view to performing the step of applying the filtering combinations; 

- the processing is spatial processing and the step of applying the 
filtering combinations comprises gain modifications applied to the coefficients of 
the initial representation; 
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- the processing is spatial and frequency processing and the step of 
applying the filtering combinations comprises the application of filters varying as a 
function of the frequency to the coefficients of the initial representation; 

- the sets of coefficients representative of the initial sound field 
5 representation and of the modified representation are sets of coefficients called 

Fourier-Bessel coefficients. 

The invention relates also to a device for processing a sound field 
representation in order to provide a representation of the sound field by applying 
this processing corresponding to at least spatial processing of the sound field, the 

10 device receiving as an input an initial representation of the sound field in the form 
of a set of coefficients which is representative of that field in time and in the three 
spatial dimensions and providing as an output a representation modified by the 
filtering combinations corresponding to the processing applied to the coefficients 
of the initial representation, characterized in that it comprises means of defining 

15 the processing, comprising: 

- means of establishing a set of directivity functions which is 
representative of the processing operation, in the form of a set of coefficients 
corresponding to the decomposition of those directivity functions on a basis of 
spherical harmonics; and 

20 - means of determining the filtering combinations on the basis of 

that decomposition, for their use by filtering means (80) receiving the coefficients 
of the initial representation as an input and providing the modified representation 
as an output. 

According to other features of the device: 

25 - the means of defining the processing operation comprise a module 

for the direct acquisition of operating data for the device in one or other of a 
plurality of formats; 

- the means of defining the processing operation comprise a module 
for the automatic interpretation of a processing instruction enabling operating 

30 data to be provided automatically in one or other of a plurality of predetermined 
formats; 

- the operating data comprise the set of directivity functions which is 
representative of the processing operation: 
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- the means of defining the processing operation comprise a 
conversion module suitable for converting the set of directivity functions from one 
or other of a plurality of predetermined formats into a set of coefficients 
corresponding to the decomposition of those directivity functions on a basis of 

5 spherical harmonics; 

- the operating data comprise a set of parameters for the processing 

operation; 

- the means of defining the processing operation comprise a module 
for calculating all or part of the set of parameters as a function of the nature of the 

1 0 processing operation: 

- the means of determining the filtering combinations comprise a 
module for calculating filtering combinations by linear combinations of the 
coefficients corresponding to the spherical harmonic-based decomposition of the 
directivity set representative of the processing operation; 

15 - the set of processing parameters comprises a set of constants and 

the module for calculation by linear combinations uses that set of constants as 
weighting coefficients; 

- the processing operation comprises at least one distortion, all or part 
of the set of directivity functions forming a set of distortion pairs corresponding to 

20 the distortion(s) to be applied to the initial representation, and the means of 
determining the filtering combinations also comprise a module for distortion by 
modification of the filtering combinations provided by the calculation module, on 
the basis of the distortion pairs, in order to provide the filtering combinations that 
are to be applied to the coefficients of the initial representation; 

25 - the processing operation comprises solely one or more distortions, all 

or part of the set of directivity functions forming a set of distortion pairs 
corresponding to the distortion(s) to be applied to the initial representation, and 
the means of determining the filtering combinations comprise solely a distortion 
module in order to provide the filtering combinations solely on the basis of the 

30 distortion pairs. 

The invention will be better understood on reading the following 
description which is given purely by way of example and with reference to the 
appended drawings in which: 

- Figure 1 is a representation of a spherical coordinate system; 
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- Figure 2 is a general flow chart of the method of the invention; 

- Figure 3 is a detailed flow chart of a first embodiment of the invention; 

- Figure 4 is a diagram of an embodiment of the step of applying the 
processing; 

5. - Figure 5 is a detailed flow chart of a second embodiment of the 

invention; and 

- Figure 6 is a block diagram of a device suitable for implementing the 
method of the invention. 

Figure 1 shows a conventional spherical coordinate system in order 
10 to clarify the coordinate system to which reference is made in the text. 

The coordinate system is an orthonormal coordinate system, having 
an origin O, and comprises three axes (OX), (OY) and (OZ). 

In this coordinate system, a position indicated x is described by 
means of its spherical coordinates (r,0,#), where r denotes the distance relative to 
15 the origin O, 0 the orientation in the vertical plane and <j> the orientation in the 
horizontal plane. 

In such a coordinate system, a sound field is known if the sound 
pressure indicated p(r,0,<j> t t) t whose Fourier transform is indicated P(r,6,<f>,j) where 
/denotes the frequency, is defined at all points at each instant t. 
20 The method of the invention is based on the use of spatio-temporal 

functions enabling any sound field to be described in time and in the three spatial 
dimensions. 

In the embodiments described, those functions are what are known as 
spherical Fourier-Bessel functions of the first kind and will be referred to 
25 hereinafter as Fourier-Bessel functions. 

In a region empty of sources and empty of obstacles, the Fourier- 
Bessel functions correspond to the solutions of the wave equation and constitute 
a basis which generates all the sound fields produced by sources located outside 
that region. 

30 Any three-dimensional sound field may therefore be expressed by a 

linear combination of the Fourier-Bessel functions, in accordance with the 
following expression: 
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In that equation, the terms Pi^ij) are defined as the Fourier-Bessel 
coefficients of the field p(r 9 0 9 <p,t), k = , c is the speed of sound in air (340ms" 1 ), 
jikr) is the spherical Bessel function of the first kind and of order / defined by 
ji{x)-^^J Mn {x) where J v 00 is the Bessel function of the first kind and of order 

v, and yr{&S is the real spherical harmonic of order / and of term m, with m 
ranging from -/ to /, defined by: 

yr(O^ = P}%os0)trg m (<f>) 

with: 

1 



trg m W = 



1 



cos( m<f>) for m > 0 
for m - 0 



Jin 

5 sin(m^) for m < 0 



10 In that equation, the P t m {x) are the associated Legendre functions 

defined by: 



with Pix) denoting the Legendre polynomials, defined by: 



15 The Fourier-Bessel coefficients are also expressed in the temporal 

domain by the coefficients pi, m (t) corresponding to the inverse temporal Fourier 
transform of the Fourier-Bessel coefficients Pi, m (j). 

In other embodiments, the sound field is decomposed on a function 
basis where each of the functions is expressed by an optionally infinite linear 
20 combination of Fourier-Bessel functions. 

Figure 2 shows a general flow chart of the method of the invention. 

In general, the method starts with a step 2 of defining a processing 
operation to be applied to a representation of a sound field and corresponding to 
at least spatial processing of that sound field. 
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Step 2 of defining the processing operation comprises a step 4 of 
establishing a set of directivity functions which is representative of the processing 
operation. 

A directivity function means a function, indicated G(O^J) t which 
5 associates a real or complex value, which is optionally a function of the 
frequency, or a range of frequencies, with each direction in space. 

Any directivity function G{0,cf>J) is decomposed into spherical harmonic 
coefficients, indicated G/, m (/), by a spherical Fourier transform. Thus, the 
directivity function G(6,<j>J) is obtained from the spherical Fourier coefficients 
10 G/,m(/) by an inverse spherical Fourier transform, in accordance with the 
expression: 

/=0/n=-/ 

In the rest of the document, the directivity functions may be denoted by 
their expression G{G^J) in direct space or in the form of spherical Fourier 
15 coefficients G h Jf). 

This set of directivity functions is provided at the end of step 4 in the 
form of a set of coefficients corresponding to the decomposition of the directivity 
functions on a basis of spherical harmonics. 

Subsequently, a step 6 of determining filtering combinations is carried 
20 out on the basis of the decomposition of the directivity functions which is carried 
out in step 4. The notion of filtering covers frequency, amplitude and phase 
filtering or a combination thereof. 

Thus, step 6 of determining filtering combinations provides filtering 
combinations corresponding to the processing to be applied in the form of filter 
25 parameters later permitting the implementation of those filtering combinations 
with the aid, for example, of a set of digital filters. 

With a view to applying the processing, the method comprises a step 8 
of establishing an initial representation of the sound field in the form of a set of 
coefficients or signals which is representative of that field in time and in the three 
30 spatial dimensions. 

In particular, this set of coefficients is a set of what are known as 
« Fourier-Bessel » coefficients as defined with reference to Figure 1. 
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The method then comprises a step 10 of applying the filtering 
combinations, defined in the definition step 2, to the coefficients of the initial 
representation as defined in step 8. 

Step 10 permits the production of a modified representation 
5 corresponding to a representation in time and in the three spatial dimensions of 
the initial sound field modified by the defined processing, in the form of a set of 
coefficients referred to as Fourier-Bessel coefficients. 

. Figure 3 shows in detail the flow chart of a first embodiment of the 
invention. 

10 In this embodiment, the processing operations applied are processing 

operations such as spatial windowing, which consists in modifying the gain of all 
or some of the virtual sources regarded as generating the initial sound field, and 
spatial convolution or spatial filtering, which consists in modifying the angular 
spread of all or some of those virtual sources, or any other analogous processing 

15 operation. 

The windowing and convolution processing may vary with the 
frequency, which enables distance and spectrum modifications to be carried out 
on all or some of the virtual sources. 

As described in Figure 2, the method starts with a step 4 of 
20 determining a set of directivity functions which is representative of the processing 
to be applied. 

In the case where the processing comprises only a single processing 
operation, such as windowing or convolution processing, a single directivity 
function, which optionally varies with the frequency, is sufficient to represent it. 

25 Advantageously, the step 4 comprises a sub-step 12 for the direct 

determination of all or part of the set of directivity functions in one or other of a 
plurality of predetermined formats. 

Additionally or alternatively, step 4 comprises a sub-step 14 for the 
input of a processing instruction which is interpreted automatically in an 

30 interpretation sub-step 16 permitting the automatic provision of all or part of the 
set of directivity functions in one or other of a plurality of predetermined formats. 
For example, the processing is used by software means in a computer and the 
processing instruction is in the form of an icon whose validation brings about the 
performance of a given processing operation by the software means. 



11 

The sub-steps 12 and 14 correspond, for example, to manual 
acquisition or also to the reception of directivity data coming from an internal 
memory or from an external device. 

As a function of the format in which the set of directivity functions is 
5 determined, step 4 may require a sub-step 18 for the conversion of the set of 
directivity functions, from a predetermined format, into a set of coefficients 
corresponding to the decomposition of the directivity functions on a basis of 
spherical harmonics. 

For each directivity function G^ m (J) of the set of directivity functions, this 
10 conversion sub-step 18 may be carried out in several ways. 

For example, the directivity function G/, m (/) is defined on the basis of 
parameters (r gy 6 q , (f> q ) specifying the position of Q virtual sources and parameters 
H q (J) specifying the frequency response of each of those sources. A radiation 
model is associated with each virtual source so that those sources as a whole 
15 are assumed to produce a sound field whose representation in Fourier-Bessel 
coefficients enables a directivity function to be associated with that field. 

For example, the radiation model associated with each source is a 
spherical wave radiation model and the directivity function G/, m (/) is obtained in 
sub-step 18 by means of the following equations: 

q = l r 9 

with 



In a variant, this calculation is carried out on the basis of parameters 
(0 g , <p q ) specifying the. directions of the Q sources. In that case, a plane wave 
25 radiation model is associated with each source and the directivity function G^Jf) 
is obtained by means of the following equation: 

In other embodiments, Gi^if) is determined on the basis of a 
combination of plane wave and spherical wave radiation models. 
30 In another case, the directivity function is defined on the basis of 

parameters representative of a parametric description of a directivity function. 
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This description is effected in direct space and associates with each direction a 
value which may be complex. Advantageously, the parameters may vary with the 
frequency. 

In that case, the coefficients of the directivity function G/^(/) are 
5 determined by using a parametric model. The parameters of the model are 
placed in a vector y(/)and the parametric model of the directivity function is 

written G(0 9 fay{f))\ The coefficients of the directivity function G/ j/n (/) are obtained 

by a spherical Fourier transform of G(0,<p,y{f)) in accordance with the equation: 

Gi, m (f) = / / G{0,4>fi(f)) yrW) sinOMd* 

J 6=0 J4>=0 

10 For example, the parametric model may correspond to a directivity 

function having a single lobe. This model is controlled by the parameters 
specifying the direction of the lobe and the parameter specifying 
the opening angle of the lobe. 

Thus, the parametric model may be a directivity function having the 
15 value 0 for any direction and having the value 1 for each direction (0,0) verifying 
the following relationship: 

(m«W)).W)>cos(S>(/)) 
where (0\(f),<f>\<f)) . (0,<p) denotes the scalar product between the directions 
(9iV),<t>m and (0,0. 

20 In that case, the coefficients Gi, m (J) of the directivity function are given 

by: 

CUP ■ 2n (J) A (J)) P '~^ A>) ~ + <»> 
In this expression P.\(x) = 1- 

In a variant, the parameters of the parametric description of the 
25 directivity function may vary with time. In that case, the parameters of the model 
are placed in a vector y(t) and the parametric model of the directivity function is 
written g(0,&Rt)). The coefficients of the directivity function gi tm (t) expressed in 
the temporal domain are then obtained by a spherical Fourier transform of 
g(0,&y(t) ) in accordance with the expression: 

».m(*)= r f 2 \(0^m)yr(0^) sine d9d^> 
30 J 0=o J(t>=o 
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The coefficients Gi^if) of the directivity function which are expressed in 
the frequency domain are then obtained by a temporal Fourier transform of the 
coefficients gi^t). 

5 In yet another case, the directivity function is defined in the form of a 

list of Q samples. The list is made up of Q pairs {((6^), v q )} where (0 q ,4> q ) is the 
direction of the sample q and where v q is the value of the directivity function in the 
direction (0 q ,$ q ). Advantageously, the list of samples may be acquired for a 
plurality of frequencies / and is then indicated {((6^), v q )}(f). 

10 For each operating frequency / the coefficients Gi im (f) of the directivity 

function are derived from the list of samples {((6^), v q )}(f). These coefficients 
are obtained by inverting the angular sampling process which permits the 
derivation of the samples from the list {((0?,^), v q )}(f) on the basis of a directivity 
function supplied in the form of spherical harmonic coefficients. This inversion 

15 may assume various forms in order to control the interpolation between the 
samples. 

In the embodiment described, the samples are placed in a vector V in 
the following manner. 

[v, v 2 —v,]' 

20 A matrix Y is constituted from the sampling directions in the following 

manner: 



^1,0,0 ^1,1 - 


1 ^1,1,0 ^1,1,1" 




-L" 




' Y \,L,L 




1 ^1,0^2,1,1" 




~L* ' 






y qs>s> y q,\- 










' Y Q,L,L_ 



where each element Y q ^ m is equal to: 
25 Y qMm =yr(0M 

It will be remembered that yr(0 q ,fa) denotes the spherical harmonic of 
order / and of term m for the sampling direction 0 q $ q considered. 

A vector G containing the coefficients G hm (f) is determined for the 
frequency /. The vector G may be determined by means of the following 
30 relationships: 
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- if the number Q of samples at frequency / is greater than the 
number of coefficients G/^(/): 

G = ( Y T Y)' 1 Y r V ; and 

- if the number Q of samples at frequency /is less than the number 
5 of coefficients G/, m (/): 

g= F T (yy T y l v. 

In the vector G, the coefficients G/, m (/) are organized in the following 

manner: 

[G 0 ,o(/) G,.-i</) G lt0 W G M (/) G^(/) - GUfi •••«]' 
10 In a variant, the list of samples may be acquired for a plurality of 

instants t and is then indicated {((0 9 ,^), v^)}(0- 

For each instant t acquired, the directivity function coefficients gi^it) 
expressed in the temporal domain are derived from the list of samples 
{((0q><f>q)> V *)X0- These coefficients are obtained by inverting the angular sampling 
15 process which permits derivation of the samples from the list {{{6 q ^ q ) y v q )}{t) on 
the basis of a directivity function supplied in the form of spherical harmonic 
coefficients. 

The directivity function coefficients G hm (f) expressed in the frequency 
domain are then obtained by a temporal Fourier transform of the coefficients 

20 gUt). 

Of course, in the case where step 12 or 16 provides the directivity 
function directly in the form of coefficients corresponding to the decomposition of 
that function on a basis of spherical harmonics, the conversion sub-step 18 does 
not modify the function. 
25 Thus, at the end of step 4, a set of directivity functions which is 

representative of the processing to be applied is provided in the form of a set of 
coefficients corresponding to the decomposition of those functions on a basis of 
spherical harmonics. 

In addition, step 2 of defining the processing operation also comprises 
30 a step 20 of establishing a set of processing parameters. 

As in the case of step 4, this set of parameters may be determined 
directly during a sub-step 22 for direct determination and/or after a sub-step 24 
for the input of a processing instruction and a sub-step 26 for the interpretation of 
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that processing instruction permitting the automatic provision of all or part of the 
set of parameters. 

The sub-steps 22 and 24 correspond, for example, to manual 
acquisition of specific parameters or also to the reception of parameter data from 
5 an internal memory or an external device. 

In particular, during step 20 of establishing a set of parameters, all or 
some of the following parameters are defined: 

- a parameter L specifying the limit order of operation of the method 
and corresponding to the desired mathematical precision limit; and 
10 - a set of constants indicated c£^ m . forming weighting coefficients 

defined as a function of the nature of the processing operation. In particular, 
these coefficients are advantageously applied to the directivity functions in order 
to obtain the above-mentioned filtering combinations. 

Some parameters may be invariant and may be directly integrated in 
15 the method. 

The parameter L enables the method to be adapted to the level of 
precision of the initial representation and, in particular, enables the method to 
which the invention relates to apply a processing operation to an initial 
representation having a high level of precision. 
20 As has been described above, the constants c\^ r ^ may be 

determined directly during step 20 by an operator, from a memory or a device 
arranged internally. For example, they are all calculated beforehand, stored in a 
memory and selected as a function of the nature of the processing operation. 

Alternatively, an instruction defining the nature of the processing 
25 operation is acquired during sub-step 24, and sub-step 26 corresponds to the 
calculation of all or part of this set of constants as a function of the nature of the 
processing operation and from relationships described later in the text. 

Of course, the constants c\$ r ^ can be determined using a 

combination of the described solutions. 
30 The parameters established during step 20, and also the coefficients 

corresponding to the directivity functions established during step 4, are used 
during step 6 of determining the filtering combinations. 
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10 



15 



20 



This determination step 6 comprises a sub-step 30 of calculating 
filtering combinations by linear combinations of the coefficients corresponding to 
the spherical harmonic-based decomposition of the set of directivity functions 
which is representative of the processing operation. 

During this first calculation sub-step 30, filtering combinations indicated 
Cf£(f) are determined on the basis of the directivity functions G lym (j) and the 

constants c^ m - which are determined as a function of the nature of the 
processing operation and which form weighting coefficients of the linear 
combinations. The filtering combinations are then determined in the following 
manner: 

r=0m"=-i 

In the above relationship, / and m denote the order and the term, 
respectively, of the coefficients of the initial representation, m' denote the order 
and the term of the coefficients of the modified representation and /", m" denote 
the order and the term of the coefficients of the directivity function. 

In the case where the processing operation is a processing operation 
of the windowing type, the constants cf;^ r ^ are determined by the relationships: 



c Z,m,J",m" 



J_ (Jm 9 (1) , J'.m f (1) \ m > 0 



J',m' (1) 



l\rn> (1) 
c i,(M",m" 



pour m = 0 



j_ / r,m' (1) _ l',m' (1) \ q 
, V2 V C /-m,/",m" C /,m,*",m">/ P OUr 



with 



Z',m' (1) 



* + Pour m« > 0 

43-5 Pour m" = 0 



with 



l',m< (2) 
c /,m,/",m" 



' C/.™' ( 3 ) + ( 3 )\ oour m' > 0 

^ l C l,m,i",m" + c l,m,l",m") P our m ^ U 

pourm' = 0 

^A^?i$-4£A>) Pourm'<0 



with 



J',m<(3) _^. m 1 V(2i + l)(2y + l)(2^+l) 
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°/-t-|m| ^r+Ki *^/"+|m"| 



2 z + z' + i" + 1 ^ c[ m| cjr' 1 cjr 1 

i /' l" rk r k-\m\ p fc / pfc'-lm'l rk <> r k"-\* 

EST V / Txfc+jfc'+Jb" M ^/ *^/' ^7" M» 

k=\m\ k'=\rn'\ k"=\m"\ 



X 



and 



and 



and 



n = 



1 pour x = y 
0 pour x 7^ y 

|m| + K| + |m"| 



n p\(n-p)\ 



This processing operation corresponds to the application of spatial 
windowing with temporal filtering since the spherical harmonic coefficients Gi >m (f) 

10 correspond to the directivity function G(0,<f>j). Likewise, the Fourier-Bessel 
coefficients P^i^f) of the initial sound field correspond to spherical harmonic 
coefficients of a directivity function P^iQ^j). Thus, the spatial windowing with 
temporal filtering results, for each frequency / in a windowing of the directivity 
function P^{04J) by the directivity function G(0,<j>j), which corresponds to a 

15 v multiplication for each direction (0,<fi) of the value of those two directivity functions, 
it being possible for the values of the two directivity functions to be complex. 

As described above, this processing applied to the sound field has the 
effect of modifying the amplitude, the distance and the spectrum of all or part of 
the set of virtual sources regarded as generating the sound field. 

20 The directivity function G(0 9 0j) may be independent of the frequency 

and is then indicated G(0,<p). The processing is then solely spatial and 
corresponds, for each frequency /, to the multiplication of the directivity function 
F*\0,<f>j) by the directivity function G(0,$. 

In the case where the processing operation is of the convolution type, 

25 the constants d^ r ^ are determined by the relationship: 
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This processing operation corresponds to the application of spatial 
filtering with temporal filtering. 

For the spherical harmonic coefficients G/, m (/) correspond to Fourier- 
Bessei coefficients of a sound field produced by a linear and angularly invariant 
sound system, in response to excitation by a reference sound field, such as a 
plane impulse wave coming from the direction (0,0). Such a response is denoted 
by the expression « spatio-temporal responses A linear sound system is 
angularly invariant when a rotation of the field which it receives induces the same 
rotation on the field which it generates. Thus, the coefficients G />m (/) correspond to 
the spatio-temporal response of such a sound system. 

The Fourier-Bessel coefficients P^i^tf) of the initial sound field are also 
regarded as spherical harmonic coefficients of a directivity function l*\0,(f>j). 
Thus, the spatial convolution with temporal filtering results, for each frequency/ 
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in a spatial convolution of the directivity function I*\0,#J) by the directivity 
function G(0,<pj) t which is regarded at each frequency / as a directional 
convolution kernel. 

In the embodiment described, the convolution is applied taking a 
5 direction (0o,<fa) into account, the convolution kernel G(0,<f>J) is then turned, 
bringing its north pole (vertical direction with 0=0 and <f> = 0) to (6b,#)). In order to 
obtain the value of the modified sound field P (r) (0,<p i /) in the direction (0',0')> the 
value in the direction (0\<f>') of the turned kernel is multiplied by the value of 
I*\0rfj) in the direction (6b,#>) and those products are added together for all of 
1 0 the directions (6b>#>). 

This processing applied to the sound field permits, for example, 
modification of the angular spread, i.e., the spatial localization precision, that is to 
say, the size of the virtual point sources regarded as generating the initial sound 
field, for each frequency f. This processing is the spatial analogue of temporal 
15 convolution by a window, the effect of which is to modify the temporal spread of a 
sound event, such as, for example, an averaging or delay effect. 

In the case where the directivity function Gi jTn {f) has cylindrical 
symmetry, that is to say, if G/^(/) = 0 for m*0, the above equations are simplified, 
and the coefficients c{£^ m .. are then determined by the relationship: 

J',m' _,l gl cm r0 / ^ 

The directivity function G{04J) may be independent of the frequency 
and is then indicated G{0,<f>). The processing is then solely spatial and 
corresponds, for each frequency / to the convolution of the directivity function 
f^\0^f) associated with the initial field by the directivity function G(0 y <p). 

25 In addition, the method of the invention also enables rotation 

processing to be carried out. In that case, however, use is not made of directivity 
functions and the method requires parameters (0,&i//) representative of a rotation 
of the sound field about the origin O. 

Thus, if the processing operation is a rotation, the filtering 

30 combinations are determined from the parameters (0,&y/) representative of the 
rotation of the sound field according to the relationship: 
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5R(rr) = partie reelle (x) 



$s(x) = partie imaginaire (x) 



In this relationship, the convention used for the angles 0, <f> and ^, 
consists in effecting a rotation through an angle y/ about the axis (OZ) then 
through an angle 6 about the axis (OY) then through an angle <p about the axis 
15 (OZ). 

Finally, the processing carried out may also correspond to neutral 
processing. In that case, the filtering combinations are all zero except for the 
values where /'= / and nC= m where they equal 1 . 

Of course, other processing operations as well as combinations of the 
20 processing operations described can be carried out thanks to the definition of the 
set of weighting constants cf;£ r%m » • 
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For example, random processing may be carried out by the random 
definition of the set of constants. 

Advantageously, the calculation sub-step 30 is carried out solely for 
the non-zero constants . 

5 The filtering combinations Cf£(f) so obtained are placed in a matrix C 

of size (L'+l) 2 by (£+l) 2 . 

In the case where the processing operation to be applied consists 
solely of rotation, convolution or filtering, the matrix C corresponds directly to the 
filtering combinations to be applied and is indicated in a general manner F so that 

10 the results of the calculation sub-step 30 form the filtering combinations that are 
to be applied to the coefficients of the initial representation. 

In the case where the processing operation corresponds to a 
combination of processing operations, such as rotations, filtering operations, 
convolutions and/or any kind of processing, the sub-step 30 is repeated a plurality 

15 of times and each matrix C is memorized in succession in order to be combined 
with the following matrices, finally to provide a matrix F comprising the 
coefficients corresponding to the filtering combinations that are to be applied to 
the coefficients of the initial representation. Thus, the matrix F is obtained by 
matrix multiplication of the successive matrices C. 

20 Each element Fj£*(f) of the matrix F represents the frequency 

response of a spatio-temporal processing filter. 

Various forms of filters and filtering methods may be used in step 10 of 
applying the filtering combinations to the initial representation of the sound field. 

In the case where the filters are parameterized directly with frequency 
25 responses, the coefficients Ff^(J) are applied directly. Advantageously, step 6 of 

determining the filtering combinations comprises a conversion sub-step 32 in 
order to determine the parameters of the filters for other filtering methods. 

For example, the filtering combinations Ff^(J) are converted into: 

- finite impulse responses f£*(t) calculated by inverse temporal Fourier 

30 transform of F^{f) , each impulse response f^\t) is sampled and then 

truncated to a length suitable for each response; or 
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- coefficients of recursive filters having infinite impulse responses calculated 
from theFfjftf) with adaptation methods. 

The filtering combinations so defined are applied during step 10 of 
applying the filtering combinations, which transforms the initial representation 
5 F^ijnif) int0 a modified representation which is indicated P^i^nif) and which 
corresponds to a modified sound field, by applying thereto the processing filters 
representative of the at least spatial processing defined in step 2. 

During step 10, coefficients /? (7 V,(0 corresponding to the coefficients of 
the modified representation are derived from the coefficients p (I) i, m (0 
10 corresponding to the coefficients of the initial representation P 40 , by the 
application of the frequency response processing filters F^{f) in the following 
manner: 




where is the Fourier transform olp ( \ m (t) and where P^i^if) is the Fourier 

15 transform of p {T) i, m (t\ these being coefficients relating to the initial representation 

and the modified representation, respectively. 

Thus, step 10 implements the conventional temporal filtering methods 

for the application of the filtering combinations provided by step 2 in a suitable 

form. 
20 For example: 

- the filtering may be effected by means of filtering methods in the 
frequency domain, such as, for example, block convolution techniques. In that 
case, the filtering combinations are supplied in the form of the frequency 
responses Ff^if). 

25 - the filtering may be effected in the temporal domain by temporal 

convolution. In that case, the filtering combinations are supplied in the form of the 
impulse responses f&Xt). 

- the filtering may be effected in the temporal domain by means of 
infinite impulse response recursive filters. In that case, the filtering combinations 

30 are supplied in the form of the recurrence relation coefficients. 

Figure 4 describes the case of finite impulse response filtering. This 
filtering initially requires the determination of a variable Tfi? corresponding to the 
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number of samples characteristic of each response f,[f(f) , which -leads to the 
following temporal convolution expression: 

/=0m=-/ r=0 

These coefficients p (7) /^(0 are representative of a sound field 
5 corresponding to the initial sound field after the application of the at least spatial 
processing which associates with a sound field produced by at least one source, 
another sound field produced by at least one source whose direction, distance 
and/or gain are thus modified. 

In another embodiment, the method of the invention also permits the 
10 application of distortion-type processing operations corresponding to a distortion 
of portions of the sound field by modifications to the characteristics of spread, 
position, amplitude of all or some of the virtual sources generating the initial 
sound field. 

This distortion processing may vary with the frequency, which enables 
15 distance and spectrum modifications to be carried out on all or some of these 
virtual sources. 

Advantageously, this processing may be combined with one of the 
processing operations as described with reference to Figure 3. 

Figure 5 describes the flow chart of the method of the invention in such 
20 an embodiment combining a first processing operation, called nominal processing 
and cprresponding to a processing operation as described above, with a spatial 
distortion processing operation which results in distortions of frequency, 
amplitude and/or phase characteristics of components of the representation of 
the initial sound field. 

25 Thus, in this embodiment, step 2 of defining the processing comprises, 

as before, the step 4 of establishing a set of directivity functions comprising the 
sub-step 12 of direct determination and/or the sub-steps 14 of acquisition and 16 
of interpretation as well as the conversion sub-step 18 enabling a set of directivity 
functions to be provided in the form of a set of coefficients corresponding to the 

30 decomposition of those functions on a basis of spherical harmonics. 

In this embodiment, part of the set of directivity functions determined in 
step 4 corresponds to N distortion pairs representative of a distortion to be 
applied. In a manner similar to the embodiment described above, these distortion 
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pairs are described in the form of a set of coefficients, indicated 
{(Hi^niftJTi^nfJ))}, corresponding to their decomposition on a basis of spherical 
harmonics. 

Each distortion pair is constituted by two directivity functions which 
5 correspond to the Fourier-Bessel coefficients of a source sound field and a target 
sound field, the distortion operation applied to the source sound field tending to 
cause it to resemble the target sound field. 

Thus, step 4 provides a set of directivity functions corresponding to the 
nominal processing and also to distortion pairs representative of the distortion to 
10 be applied in the form of a set of coefficients corresponding to the decomposition 
of the directivity functions on a basis of spherical harmonics. 

As above, step 2 also comprises the step 20 of determining the set of 
processing parameters by means of the direct determination sub-step 22 and/or 
the sub-steps 24 of acquisition and 26 of interpretation. 
15 However, in this embodiment, the step 20 of determining parameters 

also permits the determination of a parameter X which may be adjusted between 
0 and 1 and which specifies the resemblance between the final processing to be 
applied and the nominal processing determined by step 4, that is to say, the 
relative magnitude of the distortion to be effected. 
20 Thus, if X is equal to 1 , only the distortion is taken into account and if X 

is equal to 0, only the processing determined during step 4 is taken into account. 

As in the embodiment described with reference to Figure 3, step 6 of 
determining the filtering combinations comprises the calculation sub-step 30 
carried out as described above and permitting the provision of a matrix C 
25 comprising a set of filtering combinations corresponding to the application of the 
nominal processing as defined by the set of directivity functions that is 
determined in step 4. 

The set of distortion pairs is used in step 6 of determining the filtering 
combinations during a sub-step 50 of distortion by modification of the filtering 
30 combinations determined at the end of the calculation sub-step 30. 

During this sub-step 50, the filtering combinations Cf£'(f) are modified 

on the basis of the distortion pairs indicated {(Hi, m ,n(f)jr iwif))} and the 
parameter X in order to supply as an output the filtering combinations Ffi*(f) 
which are representative of the processing including the nominal processing and 
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the distortion and which are to be applied to the coefficients of the initial 
representation. 

During sub-step 50, the spherical harmonic coefficients Hi^Jf) 
representative of the N source directivity functions are placed in a matrix H of size 
5 (L+l) 2 by N in the following manner: 

H o,o, H o,o,zU) 
#,.-,.,(/) #..-., 2 (/) 
#,,o,.(/) H lfi , 2 (J) 
#,,.,.(/) #.,i. 2 (/) 

H L ,i LA (f)H L ^ 2 (f) 
H LfiA if) H L ' 0t2 (J) 
_H L ' L , x (f) H L ' La {f) 

Likewise, the spherical harmonic coefficients /T/, m ,„(/) representative of 
the N target directivity functions are placed in a matrix IT of size (L'+l) 2 by N : 

H'o,o,\(f) H'o.oaif) H'o,o,n(J) 

H\,. ul (J) i/ , ,,- 1>2 (/) n\,_ UN (f) 

H\, 0 ,i(J) H\, 0t2 (f) H\ fi , N (f) 

#'.,.,.(/) H\ A , 2 (f) H\ AA ,(f) 

0 : : : 

H'v,- l ., l (f)H' l ..- l . a U) H' L ,.,, N (f) 

H'L^if) H' L .' 9 , 2 (f) H' vfi , N (f) 

_H' LW (f) H\', L , 2 (f) H\,' L%N (f) 

The coefficients Cf£'(f) of the filtering combinations of the nominal 
processing are also placed in a matrix C of size (L'+l) 2 by (L+l) 2 : 





c 0 %°(/) cr&if) c, 0 b °c/) c,V(/). 
cd.o-'C/) ci : :i(f) c.Vo'c/) c,V(/). 
c 0 :o°C/) c.'.f.c/) c.W) c.'.fc/) • 

C^C/) C I ';J,(/) C,V 9 l (/) C };{(/) ■ 


■ C° L ; 0 . L (f)- 

■ Ci--[(f)- 

■ ci-^c/)- 


ci-m ••• 

Ci;°o(7) ••• 
Ci;K/) ••• 


ci-.i(/) 
ci;ic/) 
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••££#'(/>•• 


ct:i L \f)- 


ctX L 'if) 




Co^oo c.^c/) c.^c/) c.^c/v 

• • • « 




C £■<?(/)- 


ct;m 




.c 0 y(/) c.^Vc/) cw'(f> cff-w- 




cti L '(f)- 





Ho,o,n(/) 
#..-,.*(/) 
H X , 0 ,M 

#.,.,*(/) 
H i-l,n(J) 
H L ,o, N (f) 
H L ' L , N (f) 



The filtering combinations Fjj?(j) representative of the processing 
operation in its entirety are then determined in the course of the distortion sub- 
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step 50 corresponding to an optimization which determines a set of filtering 
combinations verifying in the best possible manner the distortion pairs and the 
resemblance between F{£(j) and C r ,£{J) . 

This optimization can be effected by means of methods such as 
5 gradient descent or matrix inversion methods and also permits the determination, 
on the virtual sources regarded as generating the initial sound field, of 
modifications that are not directly specified by the distortion pairs. 

For example, the modifications defined by the distortion pairs are 
interpolated on the portions adjacent to the directly targeted portions in order to 
10 obtain the filtering combinations, which results in progressive modifications to the 
initial representation of the sound field around the portions directly targeted by the 
distortion pairs. 

The embodiment described uses a matrix inversion method which 
leads to the determination of a matrix F during sub-step 50 according to the 
15 relationship: 

F= ( (1-A) C+ AIT H T ) ( (ia) /+ AHlf) A 
where / denotes the identity matrix of size (L+lf. 

In the case where A equals 1 and where N<(L+l) 2 , the matrix F is 

given by: 

20 F = C + (IT - C H).( H T H H T 

In the matrix F, the coefficients F,^'(f) are organized in the following 

manner: 



25 



F„V(/) F.^C/) F.yif) FiYV) 
F 0 \b l U) F'.iVC/) F.'.a-'C/) F.'.f'C/) 
F 0 \b 0 (f) Ftf(/) Ftf(f) Ftfif) 

FtMf)- F.yc/) F.'.i'c/) 



F L 0 L° L (f) 
F/A'OO 
F.'AOO 



F L °b°(J) 
FlS x {f) 



F L °;?(f) 
Fl;l x (J) 

Fl:l(f) 



Fo%\ L WF l ^WF l %- L WFA- L '(f)--- F L % L (f) ~ F L L ; 6 - L •(/)■■■ F L L 'i~ L '(J~) 
Fti°(f) F£f(f) F,f 0 ; 0 (/) F.f^C/)- F L %<J) - Ftfif) - Ftfif) 
Fo&'U) Ftf if) F^lf) F^'if)- Ftf(f):. Ftf'if)- Ftf'if)^ 



Of course, according to the interpolation method used, this 
optimization step may involve progressive or, on the other hand, sudden 
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variations in the modifications between the portions directly targeted by the 
distortion pairs and the adjacent portions. 

Subsequently, the filtering combinations Ffj?(f) are optionally 

converted during sub-step 32 and are applied during step 10 as described with 
5 reference to Figures 3 and 4. 

Advantageously, several distortions are applied in succession, the 
method then comprising a plurality of distortion sub-steps 50 which are recursive 
with respect to one another, so that each sub-step 50 receives as an input, as 
nominal processing, the filtering combinations provided by the preceding sub- 
10 step. 

In the case where the parameter X equals 0, distortion is not applied 
and the filtering combinations determined during the calculation sub-step 30 
correspond directly to the filtering combinations representative of the processing. 
Thus parameterized, this embodiment corresponds to the first embodiment 
15 described with reference to Figures 3 and 4. 

In another embodiment, only distortion is applied and the filtering 
combinations are determined solely by means of the distortion pairs in 
accordance with the relationship: 

- If AT<(Z+1) 2 : 

20 F = I+ (IT —H) ( H T H )' x H T 

- If N>(L+lf: 

F = IT H T (HH r ) 1 
In these equations, the matrices H and IT are determined as 
described above. In this case, the matrix C is not used and sub-step 30 is not 
25 carried out. Thus, step 2 comprises solely the distortion sub-step 50. 

This embodiment is a simplification of the embodiment described when 
A equals 1 and the nominal processing is neutral processing. 

It is therefore apparent that the method of the invention permits the 
application of a very wide variety of processing operations to a sound field 
30 representation having a high level of precision, owing to the use of a set of 
directivity functions which is representative of the processing operation and which 
is represented in the form of a set of coefficients corresponding to the 
decomposition of those functions on a basis of spherical harmonics. 
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In the case where the processing to be applied can be split into a 
plurality of elemental processing operations such as described above, the filtering 
combinations of the processing operation as a whole may be obtained by matrix 
multiplication of the matrices corresponding to each of the plurality of processing 
5 operations, that is to say, by arranging the filtering combinations provided by 
calculation sub-step 30 and/or distortion sub-step 50. 

This processing is thus repeated a plurality of times, the results of 
each occurrence being combined with those of the previous occurrence. 

Figure 6 shows a block diagram of a device suitable for implementing 
10 the method of the invention. 

The device 60 permits at least spatial processing of a sound field 
representation in order to provide a representation modified by the at least spatial 
processing. 

The device 60 comprises means 61 of defining the processing for the 
15 implementation of step 2 of the method of the invention. These means 61 
comprise means of establishing a set of directivity functions, which is 
representative of the processing, in the form of a set of coefficients corresponding 
to the decomposition of those directivity functions on a basis of spherical 
harmonics. 

20 Advantageously, these means comprise a module 62 for the direct 

acquisition, by an operator, from a memory or an external device, of operating 
data for the device 60, permitting, in particular, the direct determination of the set 
of directivity functions. 

For example, the acquisition module 62 is constituted by a keyboard 

25 permitting the direct acquisition of data or by an interface with an external device 
supplying directly the data necessary for the implementation of the method of the 
invention. 

The module 62 may also be used for inputting a processing instruction, 
and the means 61 then also comprise a module 64 for the automatic 
30 interpretation of that instruction, in order to provide operating data for the device 
60 and, in particular, the set of directivity functions. 

The data entered by means of the module 62 for acquisition and/or 64 
for interpretation may also have different predetermined formats and the device 
60 comprises a module 66 for converting the data between predetermined 
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formats and selected formats for the implementation of the method of the 
invention. 

In particular, the module 66 enables the set of directivity functions to 
be converted from one or other of a plurality of formats into a set of coefficients 
5 corresponding to their decomposition on a basis of spherical harmonics. 

Thus, the modules 62, 64 and 66 permit the implementation of steps 4 
and 20 of the method of the invention in order to provide a set of directivity 
functions representative of the processing operation in the form of the coefficients 
corresponding to the decomposition of those functions on a basis of spherical 
10 harmonics. 

In the case where the device 60 is used to carry out distortion 
processing, the modules 62, 64 and 66 provide directivity functions forming 
distortion pairs. 

The modules 62, 64 and 66 also permit the determination of all or part 
15 of a set of parameters necessary for the implementation of the method of the 
invention. 

In particular, those modules 62, 64 and 66 permit the definition of all or 
some of the following parameters: 

- the parameter L specifying the limit order of operation of the method 
20 and corresponding to the desired mathematical precision limit; 

- a parameter S defining the nature of the processing operation to be 

effected; 

- all or part of the set of constants indicated c££r, m « forming weighting 
coefficients and defined as a function of the nature of the processing operation; 

25 and 

- the parameter A, specifying the magnitude of a distortion to be 

effected. 

Advantageously, the device 60 comprises a memory enabling all or 
some of those parameters and especially the set of constants to be stored. For 
30 example, such a memory comprises a plurality of sets of constants, each 
corresponding to a processing operation, the set of constants used being 
selected from that plurality as a function of the parameter S representative of the 
processing operation. 
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Alternatively, the means 61 of defining the processing operation 
comprise a calculation module suitable for determining all or part of the set of 
constants as a function of the parameter S representative of the nature of the 
processing operation. 

5 All of these operating data are introduced into means 70 of 

determining filtering combinations, which means 70 are included in the means 61 
of defining the processing operation and are suitable for implementing step 6 of 
the method. The means 70 comprise a module 72 for calculation by the 
application of linear combinations permitting the implementation of the calculation 

10 sub-step 30 of the method. The module 72 receives the set of directivity functions 
as an input in order to provide filtering combinations by linear combinations of the 
coefficients corresponding to the spherical harmonic-based decomposition of the 
sets of directivity functions representative of the processing operation. 

Advantageously, the constants provided by the modules 62, 64 and 66 

15 form weighting coefficients and are used as coefficients of the linear 
combinations in the module 72. 

In the case where the device 60 also uses distortion processing, it also 
comprises a module 74 for distortion by modification of the filtering combinations 
provided by the module 72, implementing sub-step 50 of the method. This 

20 module receives as an input the distortion pairs representative of the distortion to 
be applied, as well as the parameter X representative of the distortion effort used 
as described above. 

In the case where only distortion is applied, the module 72 is not 
involved in the use of the device and only the distortion module 74 is necessary. 

25 Thus, the means 70 provide a set of coefficients corresponding to 

filtering combinations representative of the processing to be applied. 

Optionally, these filtering combinations are introduced into a 
conversion module 76 suitable for providing filter parameters from the filtering 
combinations by implementing sub-step 32 of the method as described above. 

30 These filter parameters or, directly, the filtering combinations are 

applied to filtering means 80 which implement step 10 of the method and receive 
as an input the initial representation of the sound field in the form of a set of 
coefficients which is representative of that field in time and in the three spatial 
dimensions in order to apply the filters determined in order to provide a modified 
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representation corresponding to the initial sound field modified by the processing 
defined above. 

For example, the filtering means 80 are digital filtering means applied 
by a signal-processing processor to a digital representation of the initial sound 
5 field. 

In other embodiments, these filtering means are gain-controlled 
analogue filtering means applied to an analogue representation of the sound field 
in the form of a set of signals each corresponding to a coefficient of the initial 
representation of the sound field. 

10 For example, the device 60 may be in the form of a computer equipped 

with software means for calculation and hardware means for acquisition and 
reproduction, in order to permit the use of the data-input means 62 and the 
acquisition of the initial representation of the sound field and the output of the 
modified representation of the sound field, the various elements of the device 60 

15 being formed by combinations of software means and hardware means. 

In other embodiments the device 60 corresponds to a specific 
acquisition card dedicated solely to the performance of sound processing or is in 
the form of a component programmed for the implementation of those functions 
only. 

20 Advantageously, the device 60 is-., integrated in a sound system 

comprising means of acquiring a sound field and means of representing that 
sound field in the form of a set of coefficients in time and in the three spatial 
dimensions in order to provide the initial representation of the sound field. Such a 
system also comprises a reproduction unit controlled by control means receiving 

25 as an input the modified representation of the sound field and controlling a 
reproduction unit in order to obtain the reproduction of the modified sound field. 



